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ABSTRACT 

I examine the hypothesis that many of the high velocity white dwarfs observed 
by Oppenheimer et al. (2001) are the remnants of donor stars from binaries that 
produced type la supernovae via the 'single degenerate' channel. If this channel 
is a significant contributor to the Galactic SN la supernova rate, then the local 
density of such remnants with V± > lOOkm.s -1 could be as high as 2 x 10~ 4 pc~ 3 , 
comparable to the densities found by Oppenheimer et al. This white dwarf 
population differs from others in that it is composed exclusively of single stars. 

Subject headings: binaries: close — stellar dynamics — supernovae: general - 
Galaxy: kinematics and dynamics — white dwarfs 

1. Introduction 

The notion that white dwarfs may play a part in the Galactic dark matter has been 
actively investigated in the last few years, motivated primarily by the microlensing obser- 
vations (Alcock et al. 1997, 2000). Recently there have been claims that this population 
has been detected directly (Oppenheimer et al. 2001). However, neither of these claims is 
without detractors. The microlensing optical depth to the LMC may also receive a signifi- 
cant contribution from LMC-associated populations as well (Sahu 1994; Zhao & Evans 2000; 
Alcock et al. 2001), while the high velocity white dwarf sample may possess a thick disk 
component (Reid, Sahu & Hawley 2001; Reyle, Robin & Creze 2001; Koopmans & Blandford 
2001; Hansen 2001). 

In this paper I wish to investigate a potential source of high velocity white dwarfs 
different from the usual population; one which may appear in the searches for halo white 
dwarfs. The proposed population results from the formation of white dwarfs from stars in 
the range 1.3 — 3 M Q which act as the donor stars in close mass-transfer binaries in which 
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the accretor is a white dwarf. Observationally, these binaries are thought to be the origin of 
the Super-Soft X-Ray sources. These binaries may also contribute to the type la supernova 
rate in spiral galaxies. If so, the explosion of the accretor will release the donor star with its 
pre-supernova orbital velocity. Thus, the white dwarf remnants of these donors (hereafter 
DR, or 'donor remnants') will be a population of high velocity white dwarfs. It is the origin 
and implications of this population that we investigate below. 



2. The Single Degenerate Supernova Scenario 

The working hypothesis of this paper is that the observed Type la supernovae are the 
result of the so-called 'single degenerate' scenario. In this instance, the pre-supernova white 
dwarf accretes material from a close, non-degenerate binary companion. This is only one of 
the currently viable models for Type la supernovae (e.g., Branch et al. 1995). Indeed, it 
cannot be the only one as it is restricted to younger stellar populations (see below). However, 
if most close Super-Soft X-Ray sources lead to Type la supernovae, they will produce a rate 
comparable to the estimated Galactic rate (Kahabka & van den Heuvel 1997; Kahabka 1999). 

In order for the white dwarf to grow the accretion rate must meet several severe con- 
straints (Rappaport et al. 1994; Kahabka & van den Heuvel 1997). First, the accretion rate 
must be high enough so that steady burning is possible, otherwise nova eruptions remove 
the accumulated mass and the white dwarf never grows (and may even experience a net 
loss of mass). On the other hand, the accretion rate must be low enough that the accretor 
envelope does not swell to the point that it engulfs the companion, which would result in a 
common-envelope episode, effectively ending the mass transfer. 

Following the Rappaport et al. (1994) model, we find that the accretion rate is limited 

by 

M < 8 x lO^Meyr" 1 (M wd - O.5M ) (1) 
M > 1.3 x l(T 7 M M3 d 6 (2) 

For a donor that is only slightly evolved, so that the envelope still has a significant 
radiative component, the mass transfer occurs on the thermal time of the envelope (if the 
donor mass is larger than the accretor mass). Thus, we approximate the accretion rate as 

M = 0.5 [M d (0) - M wd (0)] /r KH (3) 

where tkh = 3.1 x 10 7 yrsM d " 2 is the Kelvin-Helmholtz time. Thus, to meet the requirements 
for stable mass transfer, we require 1.3M < M d < 3.3M . When the type la supernova 
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occurs, the remnant of the donor will emerge with the pre-supernova orbital velocity. The 
fact that M d > 1.3M means that all the kicked stars end up as white dwarfs. Even after 
taking into account the significant mass loss, Langer et al. (2000) find that all the donor 
remnants have M > 1M Q (see their Figure 18). Thus, this scenario produces a high velocity 
population composed exclusively of white dwarfs. 

The velocity distribution depends on the separation of the binaries when the explosion 
occurs. The population synthesis of Super-Soft Sources by Rappaport et al. (1994) reveals 
binaries with orbital periods in the range 0.3 - 2 days. This means orbital velocities of 
~ 100 — 200 km.s^ 1 , depending on the mass ratio of donor and white dwarf at ignition. This 
is lower than the velocities found by Canal, Mendez & Ruiz-Lapuente (2001) because this is 
the absolute space velocity of the donor, rather than the velocity relative to the supernova 
remnant centre (a factor M wd /(M wd + M donor ) smaller). 

3. The Galactic Population 

To determine the importance of this DR population to high velocity white dwarf studies, 
we need to know both what the overall birthrate is and also how the kick velocities modify 
the original thin disk kinematics. 

Estimates of the local type la supernova rate (Tammann, Loffler & Schroder 1994; 
van den Bergh & McClure 1994) yield ~ 4 x 10~ 3 yr _1 for the Galaxy, suggesting a total 
population of ~ 4 x 10 7 remnants in the Galaxy. This is an estimate of the current rate. 
One may also estimate a rate based on the total amount of Iron in the Galaxy (~ 1O 8 M 
assuming solar proportions and 9 x 10 10 M Q in the Galactic disk). It depends on what one 
assumes for the average output per supernova and the relative contribution of core collapse 
and thermonuclear supernovae. Low estimates (~ 9 x 10 7 type la) are obtained using the 
results of Arnett, Schramm & Truran (1989) (~ 0.7M Q per event and ~ 50% of the global 
Iron from type II supernovae) while high estimates (~ 2.8 x 10 8 type la) are obtained using 
the numbers from Qian & Wasserburg (2002) (~ O.3M per event and ~ 1/3 from type II). 
So, we will consider a range for the global numbers of remnants born ranging from 4 x 10 7 
to 2.8 x 10 8 . 

There is also no quarantee that every supernova la arises from the kinds of progenitor 
system assumed here. Branch et al. (1995) describe the Super-Soft sources as a promising 
candidate for the progenitor populations in young stellar systems such as spiral galaxies, 
but not for elliptical galaxies (in which the > 1.3M objects should have evolved to white 
dwarfs by now). Hillebrandt & Niemeyer (2000) favour the single degenerate scenario above 
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the double degenerate or merging white dwarf scenario. Whether a realistic birthrate can 
be achieved depends on the uncertain details of population synthesis studies (Yungelson et 
al. 1996; Di Stefano & Nelson 1996; Kahabka & van den Heuvel 1997). I will subsume this 
uncertainty into the general rate uncertainty described above. 

To describe the effect of the kick velocity, I simulate the DR population by applying 
randomly oriented kicks to a thin disk progenitor system. I use the simple Galactic model 
of Paczynski (1990), with bulge, halo and disk potentials 



where j = b, d,M b = 1.12 x 1O 1O M , M d = 8.78 x 1O 1O M , M h = 5 x 1O 1O M , a h = 0,a d = 
4.2 kpc, b h = 0.277 kpc, b d = 0.198 kpc and r h = 6 kpc. I assume stars are born from 
exponential density distributions in both R and z, with scale lengths, 2.8 kpc and 0.35 kpc 
respectively. I assume the kick velocities are evenly distributed between lOOkm.s -1 and 
200km.s _1 , and randomly oriented. The epoch of the supernova event is distributed evenly 
over the last 10 Gyr. The resulting DR proper motion distribution is shown in Figure 1. 

The result of these simulations is a total density near the sun 



where N tot is the total number of type la supernovae that have exploded in the Galaxy over 
the last Hubble time. To compare directly with the estimated local density from Oppen- 
heimer et al. (2001), we must take only that fraction (~ 0.34) which passes the proper motion 
cut Vj_ > lOOkm.s -1 in the solar neighbourhood. Thus, the range of possible densities in 
the Oppenheimer sample, depending on the global supernova rate, is n ~ 2.6 x 10~ 5 pc~ 3 - 
1.8 x 10 -4 pc~ 3 . Thus, at the very least the DR density is comparable to the expected 
spheroid contribution 1 and the value at the upper end of the allowed range would represent 
a major contributor to the density found by Oppenheimer et al. 



1 This can be easily understood as the estimated total stellar mass in a stellar population per SN la is 
~ 7OOM (Ruiz-Lapuente & Canal 2001), similar to the amount of thin disk stellar mass per spheroid star 
in the solar neighbourhood. 







(6) 



- 5 - 



3.1. Telling the difference 

The origin of the Oppenheimer et al. (2001) sample has been a subject of some dispute. 
How do we determine whether the scenario described above is important or not? The 
individual properties of the DR white dwarfs themselves do not reflect much about their 
origins. However, we may potentially test this scenario by considering the unique aspects of 
the population as a whole. The most important is that the DR population is composed solely 
of single stars. The fact that the kick velocity is actually the orbital velocity of a disrupted 
close binary excludes the possibility that any of the DR have a binary companion. Close 
binary companions are excluded by the fact that the progenitor system is a mass-transferring 
close binary, i.e., there is nowhere one could place a third body that would be dynamically 
stable and furthermore not interfere with the mass transfer process. Similarly, a distant 
third body bound to the original binary would no longer be bound after the supernova. This 
follows from the fact that the third body must have a separation at least twice that of the 
close binary (to be stable) and a consequently smaller orbital velocity. 

The binary frequency amongst the local white dwarfs is 25% (Holberg, Oswalt & Sion 
2002). One may also estimate this by the Vmax-weighted binary contribution to the proper 
motion-selected sample of Liebert et al. (1988). Using the weights of Leggett, Ruiz & Berg- 
eron (1998), we find 8% of the white dwarf number density lies in close double degenerate 
stars, whose absolute magnitudes are too bright to be from a single C/O white dwarf. An- 
other 15% is contributed by white dwarfs in wide, common-proper motion binaries with 
main sequence stars brighter than the faint end of the white dwarf sample. The fraction 
may be higher if there are a lot of faint, low-mass companions, but the estimate of 23-25% 
is a conservative limit. Thus, the acquisition of accurate trigonometric parallaxes for the 
Oppenheimer sample should allow us to determine the fraction with anomalously bright lu- 
minosities for their colours, i.e., unresolved double degenerate binaries 2 . Further searches 
for common proper motion companions (as bright or brighter than the white dwarfs them- 
selves) will constrain the wide binary fraction. With sample sizes of 100 or more objects, 
this hypothesis is eminently testable with current capabilities 3 . 

However, the situation can be complicated if the DR population is mixed with a compa- 
rable density of spheroid white dwarfs or even thick disk white dwarfs (Reid, Sahu & Hawley 



2 In principle, single Helium core white dwarfs could also show this signature. However, Langer et al. 
(2000) find that all supernova binaries leave donors with M > 1M , which should produce C/O core white 
dwarfs. 

3 We stress, however, that the binary fraction quoted refers to the I4iax-weightcd density, not the simple 
number count. 
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2001; Reyle, Robin & Creze 2001). Estimates of the densities from these populations are 
somewhat uncertain, but we can hope to seperate them by examining the age distributions. 
In particular, both these populations are thought to be the result of old bursts of star forma- 
tion early in the history of the Galaxy. On the other hand, the DR population is drawn from 
a stellar population that is continuously forming stars. Thus, even if all the other physics 
is the same, the different star formation histories lead to different white dwarf luminosity 
functions. 

Consider now a white dwarf population that contains a total density composed equally 
of thick disk or spheroid remnants (i.e. a burst population) and DR (a continuously produced 
population). Figure 2 shows the fraction of the white dwarfs in each luminosity bin that 
belong to the burst population, i.e. the fainter bins are dominated by the burst dwarfs and 
the brighter bins by the DR. This demonstrates how we might disentangle the contributions 
from the two populations. If the burst population has a binary fraction of 25% and the 
DR population has no binaries, then the binary frequency of the combined population, for 
My < 14, is ~ 1/3 x 0.25 ~ 8%. Thus, even if the total contributions are similar, the binary 
fraction amongst the bright white dwarfs is measurably reduced. 

In fact, the Oppenheimer et al. (2001) sample contains a lot of hot white dwarfs, 
enough so that Hansen (2001) has argued that the luminosity function contains a significant 
fraction of stars with a star formation history more consistent with the thin disk than an 
old burst population. If this is true, then the binary fraction of this sample should be low. 
We note that two other scenarios have been suggested to address this problem. Koopmans 
& Blandford (2002) and Davies, King & Ritter (2001) have also suggested mechanisms by 
which stars can be ejected from the thin disk to form an extended population. However, 
the one described here is the only one guaranteed to consist entirely of white dwarfs. As 
such, a control study of binary fractions in another stellar type (e.g. M dwarfs) with similar 
kinematics, can potentially also be used to test this scenario. 

The asymmetric drift of a population can provide some information about the origins. 
Using the simulations we compare the mean tangential velocity at the solar neighbourhood 
to the circular velocity, and we find an asymmetric drift V a ~ V ciic — < V<p >~ — SSkm.s^ 1 . 
So the post-kick kinematics of this population are closer to that of a thick disk population 
than a halo population. Similar drifts are produced by the Koopmans & Blandford (2002) 
& Davies et al. (2001) scenarios. 
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3.2. Microlensing Optical Depth 

It is of interest that the DR population is composed of single stars only. As such, it 
circumvents a very popular argument that ascribes the dominant microlensing contribution 
to lens populations associated with the Magellanic clouds themselves. The argument in 
favour of this is that the binary lenses so far detected (which provide extra information - 
namely the lens proper motion) are located in the clouds (Albrow et al. 1999; Afonso et al. 
2000, Alcock et al. 2001). In the case of the DR, no such events are even possible. 

One can estimate the microlensing optical depth from our simulations 

[ D -4wG D d (D, - D d ) Jn 

where D s = 50 kpc. From the simulations, the integral ~ 0.0017, so that r = 6.9 x 10~ 10 , 
i.e. a negligible contribution to the observed lensing optical depth 1.2 x 10~ 7 . The essential 
problem is just that the total mass in the population is very small. A similar estimate for 
the DR population of the LMC itself yields ~ 2.1 x 10~ 9 . So, this population, either in the 
Galaxy or the LMC, is unlikely to be a viable solution to the microlensing problem. 



4. Conclusion 

If the majority of Type la supernovae in the Galaxy are the result of single degenerate 
binaries with slightly evolved main sequence donors, the post supernova kinematics of the 
donors leads to a population of high velocity white dwarfs. These can potentially explain 
the interesting observations of Oppenheimer et al. (2001). Not only is the density sufficient 
to explain the observations, but the age distribution is more characteristic of the thin disk 
than the thick disk, as is suggested by the Oppenheimer et al. luminosity function (Hansen 
2001). 

The observational test of this proposal is to search for evidence of binarity amongst the 
high velocity white dwarfs, by obtaining accurate trigonometric parallaxes and searching for 
common proper motion main sequence companions, particularly amongst the hotter white 
dwarfs, where the contribution of the spheroid or thick disk is smaller. Although the DR 
represent a new population of high velocity white dwarfs and might be considered interesting 
from the point of view of microlensing, the total mass involved is quite small and the expected 
contribution to the microlensing optical depth from this population is negligible unless the 
Galaxy retains only a small fraction of the total iron produced by Type la supernovae. 
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Fig. 1. — The distribution of proper motions for donor remnants within 1 kpc of the sun, 
based on the simulation. We find that 34% of the systems lie above the 100km. s" 1 cut 
applied by Oppenheimer et al. 

Fig. 2. — The fraction of stars Af burst /(iV burst + iV DR ), in any given My bin, drawn from the 
burst population. The total integrated number of stars in both populations is the same and 
the age of the burst is assumed to be 12 Gyr. The DR population has been forming stars at 
a constant rate over the same period. Thus, we see that the burst stars pile up at the faint 
end and are less common at the bright end. 
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